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Abstract 

We introduce isospin dependence in the cluster recognition algorithms used in the Quantum 
Molecular Dynamics model to describe fragment formation in heavy ion collisions. This change 
reduces the yields of emitted nucleons and enhances the yields of fragments, especially heavier 
fragments. The enhancement of neutron-rich lighter fragments mainly occurs at mid-rapidity. 
Consequently, isospin dependent observables, such as isotope distributions, yield ratios of n/p, 
t/ 3 He, and isoscaling parameters are affected. We also investigate how equilibration in heavy ion 
collisions is affected by this change. 
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The productions of nucleons, light particles and intermediate mass fragments (IMF) in 
heavy ion collisions(HICs) provide unique experimental information about the equation of 
state(EoS) as well as other equilibrium and non-equilibrium properties of nuclear mat- 
ter far from its normal state 5Hlll|. Interpretation of the experimental results requires 
comparisons of data to predictions from transport models which take into account reaction 
dynamics. Transport models track the time evolution of a nuclear reaction, and separately 
treat the nuclear EoS through the mean field part and the nucleon-nucleon (NN) collision 
cross sections through collision part. The influence of the mean field, both isospin depen- 
dent and independent parts, and the NN collision cross sections has been studied using 



both the Bo 
approaches 



tzmann-Uehling-Uhlenbeck (BUU) and Quantum Molecular Dynamics (QMD) 



12H14|. These studies a 
description of reaction dynamics (2), 



so point out the importance of cluster formation in the 



13 



16]. 



In the QMD model, fragments are formed due to A-body correlations caused by the 
overlapping wave packets and are identified by the cluster recognition method which also 
plays important roles on the final observables. At any time of the reaction process, frag- 
ments can be recognized by a minimum spanning tree (MST) algorithm^, In this 
algorithm, nucleons which have a neighbor within a distance of coordinate and momentum 
of \ri — rj\ < Rq and \pi — Pj\ < Pq belong to a fragment. Here, the and pi are the cen- 
troids of the wavepacket for i th nucleon in their spatial and momentum space. Rq and Pq are 
phenomenological parameters determined by fitting the global experimental data, such as 
the intermediate mass fragments multiplicities 3, [l3 • They should roughly be in the range 
of nucleon-nucleon interaction, 
are about 3.5 fm and 250 MeV/c 



ypical values of R and P used in the QMD approaches 
3, Q • This approach has been quite successful in explain- 



ing certain fragmentation observables, such as charge distributions of the emitted particles, 

dnn 

intermediate mass fragments multiplicities [3, I17H19]. yield ratios of free neutron to proton 
(n/p), and the double n/p ratios in heavy ion collisions 15]. On the other hand, the MST 
method fails to describe other details in the production of nucleons and light charged parti- 
cles 0, 3, 



20j. For example, while the yields of Z=l particles are overestimated, the yields of 



Z=2 particles are underestimated partly due to the strong binding of a particles. Enhance- 
ments of the production of neutron- rich isotopes observed in isoscaling 21], dynamically 



emitted heavy fragments 
(LCP) at mid-rapidity 



22] in neutron-rich HICs, and neutron-rich light charged particles 



23] have not been very described well by simulations using transport 



models. Furthermore, most transport models predict more transparency than that observed 



experimentally in central collisions at intermediate energy [12J, |24| due to insufficient pro- 
duction of fragments in the mid-rapidity region. Previous studies show that these problems 
can not be resolved by only changing the mean field or nucleon-nucleon cross section in 
transport models. 

There have been many attempts to improve the MST algorithm. More sophisticated algo- 



rithms such as the early cluster recognition algorithm (ECRA) 25j, the simulated annealing 



clusterization algorithm (SACA) 26|, 27], and the minimum spanning tree procedure with 



binding energy of fragments (MSTB) 28j have been developed to provide better description 
of the IMF multiplicities or the average Z max of the fragments. However these algorithms 
do not address the lack of isospin dependence in cluster recognition which is the main focus 
of this paper. 

In the regular MST method, neutrons and protons are treated equally in the cluster 
recognition process, a constant maximum distance for the nucleons to be defined as belonging 
to a cluster is Rq" 1 = Rq P = Rq P = R ~ 3.5/m. In order to investigate the effects 
of introducing the isospin dependence in the cluster recognition, the iso-MST, we choose 
different values of R% n , R% p , and R P P . We found that R% n = Rq P = Qfm and R P P = 3fm 
give the largest effect on the heavy ion collisions observables in this study. These values are 
consistent with the Ro values ranging from 3 to 6 fm used in the QMD simulations in order 
to reproduce the cluster distribution reasonably. The larger distances of Rq u and Rq P take 
into consideration of the properties of neutron-rich nuclei, such as neutron skin or neutron 
halo effect. Due to the long range repulsive Coulomb force between nucleon in the cluster, 
a smaller R^ value is adopted. P Q = 250MeV/ c is the same in both the MST and the iso- 
MST methods because nucleons with large relative momentum are no longer close together 
in coordinate space after some time. In this paper, we compare the results calculated with 
the MST and the iso-MST methods to investigate effects of introducing isospin dependence 
in the identification of cluster formation on isospin sensitive observables. We also investigate 
how these cluster recognition methods affect the observable on system equilibrium. All of 



our studies are based on the ImQMD05 code 
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19]. 



Within the ImQMD05 approach, nucleons are represented by Gaussian wavepackets and 
the mean field acting on these wavepackets is derived from an energy functional with the 
potential energy U that includes the full Skyrme potential energy without the spin-orbit 



term: 

U = U p + Umd + Ucoul (1) 

and U cou i is the Coulomb energy. The nuclear contributions are represented in a local form 
with 

U P:m d = J u P:md d 3 r (2) 

and 

2 Po »7 + 1 Po 2 Po 

. 9sur,iso rv7 / \ -i 2 

+ [V(p„ -pp)] 

Po 

+y(f)^ 2 P + ^^ (3) 

1 P° Po 

Here, 5 is the isospin asymmetry. 5 = (p n — p p )/(p n + P P ), Pn and p p are the neutron and 
proton densities, respectively. A symmetry potential energy density of the form ^ (p^) 7 ^ 2 P 
is used in the following calculations. In the present work, we adopt the asy-soft value 



m 



291 ] . The energy density 



of 7j = 0.5 which provides the better agreement with data 
associated with the mean-field momentum dependence is represented by 

1 1 r 

Umd = / d 3 p 1 d 3 p 2 f N1 (f 1 )f N 2(p2) 

2p l07T D J 

1.57[ln(l + 5 x 10" 4 (Ap) 2 )] 2 (4) 

where /jv are nucleon Winger functions, Ap = \p[ — p~2.\i the energy is in MeV and momen- 
tum is in MeV/c. The resulting interaction between wavepackets is described in Ref.[?j]. 
In this work, the a = -356MeV, f3 = 303MeV, r\ = 7/6, and g sur = 19A7MeVfm 2 , 
g SU r,iso = —11.35MeVfm 2 , C s = 3B.19MeV, and g pT = OMeV. These calculations use 
isospin-dependent in-medium nucleon-nucleon scattering cross sections in the collision term 



17H19|. 



and the Pauli blocking effects are described in 

Fig.l (a) shows results of charge distributions from the central collisions of 124 Sn + 124 Sn 
at b=2fm at 50 MeV per nucleon beam energy calculated in the MST (open circles) and 
the iso-MST (solid circles), respectively. In the current paper, we adopt the convention 
that open circles are results obtained with the MST method, and solid circles represent 
results using the iso-MST method except in Figure 4. Fig. 1(b) plots the yield ratios from 
the iso-MST and the MST algorithms, i.e., Y(iso — MST) /Y (MST). Since larger values 
of Rq u and FQ P are adopted, more nucleons, especially neutrons, are included in clusters. 
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FIG. 1. (Color online) (a) Charge distribution of 124 Sn + 124 Sn at beam energy of E/A = 50MeV 
for b=2fm for ji = 0.5. The solid circles are results from the iso-MST algorithm and the the open 
circles for the same reaction system. (b)Y(iso — M ST) /Y (M ST) ratio obtained from the same 
reaction system. 



As a consequence, the multiplicities of Z=l particles are reduced by about 16%, and the 
production of heavier fragments with Z > 12 are strongly enhanced. This may explain the 
strong enhancement of dynamically emitted heavy fragments observed in neutron-rich heavy 
ion collisions 22J. Even for intermediate mass fragments with Z=2-12, the multiplicities are 
larger by ~3%. In all cases, the enhancement is larger in neutron- rich heavy ion collisions. 

In order to understand the enhancement of neutron-rich particles, we analyze the rapidity 
distributions of the particles emitted in central collisions of 124 S'?7.+ 124 S'n reactions at b=2fm 
and incident energy of E/A = 50MeV using both the MST (open circles) and the iso-MST 
(solid circles) algorithms. In Fig.2, we present the rapidity distributions of n, p,and their 
ratios in (a), (b) and (c), respectively. The calculations with the iso-MST method (solid 
circles) reduce the yields of both neutrons and protons over all rapidity regions relative to the 
results obtained with MST case, especially at mid-rapidity region. However, the yield ratios, 
Y(n)/Y(p) (Figure 2 (c)) are enhanced at mid-rapidity and become smaller at forward and 
backward rapidity regions in the iso-MST case than that in the MST case. 

In Figure 2 (d), (e) and (f), we plot the rapidity distributions of t, 3 He and their ratios, 
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FIG. 2. (Color online) The rapidity distribution of (a) the neutron yields, (b) the proton yields 
(c) the yield ratios, Y(n)/Y(p). (d) the t yields, (e) the 3 He yields and (f) the yield ratios, 
Y(t)/Y( 3 He). 



Y(t)/Y( 3 He), respectively. The t and 3 He yields are much lower than that of n and p. 
Note that the scales of the y-axis for Fig. 2 (d) and (e) are a factor of 10 smaller than 
those in Figure 2(a) and (b). Unlike the neutron yields, the triton yields (Fig. 2d) obtained 
with the iso-MST method are enhanced at mid- rapidity, and remain nearly the same in 
the projectile or target regions relative to the results obtained with MST method. On the 
other hand, there is very little difference in the yields of proton-rich 3 He (Fig. 2e) obtained 
with the iso-MST and the MST methods over the whole rapidity regions. As a consequence, 
the yields ratios, Y(t)/Y( 3 He), are much enhanced as shown in Figure 2 (f). The values 
of Y(i)/Y( 3 He) obtained with the iso-MST method increase to about 4.0 at mid-rapidity, 
much higher than the Y(n)/Y(p) ratios. Two factors contribute to the larger Y(t)/Y( 3 He) 
ratios at mid-rapidity. First, the emitted nucleons at mid- rapidity have lower kinetic energy, 
because most of them pass through the expansion phase and dissipate their kinetic energy. 
At the end of the simulations, the emitted nucleons with lower kinetic energy are closer to 
the other clusters. Thus, the iso-MST algorithm with larger Rq 11 and Rq P values will reduce 
the yields of nucleons at mid-rapidity as more nucleons close to the clusters are absorbed into 
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fragments. Furthermore, production of neutron-rich fragments is enhanced at mid-rapidity 
due to larger R^ l l np values adopted in the iso-MST case. The calculated results in Fig. 2 
clearly demonstrate that cluster recognition method changes the yields of nucleons and light 
charged particles, and neutron-rich particles emitted at mid-rapidity in heavy ion collisions. 

We also analyze the influence of the isospin dependent cluster recognition method on 
isospin dependent observables, such as R(n/p) = Y(n)/Y(p), DR(n/p), R(t/ 3 He) and 
DR{t/ 3 He) at lower and higher kinetic energy regions. In Table I, we list the yield ratios of 
n/p and t/ 3 He with Ek < 40MeV and Ek > 40MeV for the neutron rich reaction system 
124 Sn+ 124: Sn. It is clear that the iso-MST case increases the values of R(n/p) at lower kinetic 
energy relative to the MST case. To reduce systematic uncertainties in experiments, we also 
construct the double ratios DR(n/p) obtained from two reaction systems 12A Sn ± 124 Sn, 
U2 Sn + 112 Sn at b=2fm using both the iso-MST and the MST algorithms. Here 

= [y(n)/y( P )], MS „ +ms „ 

[Y(n)/Y(p)\ii2 Sn+ ii2 Sn 

For nucleons with low kinetic energy (< AOMeV), DR(n/p) obtained in the iso-MST case 
is larger by about 8.3%. For nucleons with high kinetic energy (> AOMeV), the difference 
between the DR(n/p) values obtained from both the MST and the iso-MST cases is much 
less. As most of the emitted triton and 3 He have kinetic energies less than 40 MeV per 

TABLE I. The calculated results of R(n/p), R(t/ 3 He) obtained for 124 5n + 124 Sn at b=2fm, the 
results of DR(n/p), DR(t/ 3 He) and a are obtained with an angular gate of 70° < 9 cm < 110° from 
two reaction systems 124 Sn+ 124 Sn, 112 Sn+ 112 Sn at b=2fm. 

Observable^ MST iso-MST 

R(n/p) Ek <40Mev( 12i Sn) 2.131 ± 0.005 2.302 ± 0.006 

R{n/p) Ek>A0MeV { l24 Sn) 1.044 ±0.011 1.041 ±0.011 

DR{n/p) Ek<mMeV 1.514 ±0.005 1.641 ± 0.006 

DR{n/p) Ek>mMeV 1.803 ±0.029 1.852 ±0.030 

R{t/ 3 He)i24 Sn 3.031 ±0.026 3.497 ± 0.028 

DR(t/ 3 He) Ek/A<40MeV 1.57 ±0.02 1.64 ±0.02 

Isoscaling parameter a 0.19 0.22 

a All the results are obtained with ji = 0.5 case. 
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nucleoli, we list the R(t/ 3 He) and DR(t/ 3 He) values with E k /A < AOMeV in the fifth and 
sixth row of Table I. The iso-MST algorithm enhances the R(t/ 3 He) and DR(t/ 3 He) values 
more than R(n/p) and DR(n/p). 
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FIG. 3. (Color online) Calculated primary isotope yields of Li, Be, B, C, N and for 124 S , n+ 124 Sn 
with an angular gate of 70° < 9 c . m . < 110° for 7, = 0.5 at b=2fm. 

With larger R™ and Rq P values adopted in the iso-MST method, more neutron-rich 
isotopes are produced. This effect is further enhanced in the neutron-rich reaction system. 
Fig.3 shows the isotope distributions of primary fragments with Z=3 to 8 predicted by the 
ImQMD05 calculations over the angle region, 70° < 9 < 110° for 124 Sn + 124 Sn system. 
The solid circles are results obtained with the iso-MST method, and open circles are results 
obtained with the MST method. 

The isoscaling relationship, which is constructed using isotope yields Yi(N, Z) with neu- 
tron number N and proton number Z from two different reaction systems denoted by the 
index i, obeys a simple relationship 

R 21 (N, Z) = Y 2 {N, Z)/Y 1 (N, Z) = Cexp(aN + f3Z) (6) 

Here, C is an overall normalized factor and a and f3 are isoscaling parameters. The isoscaling 
parameter a can be used to estimate the enhancement of neutron rich isotopes, and has been 

8 



used as a probe to study the density dependence of symmetry energy |2l|, [30j, |3l|, |34| • We 
analyze the isoscaling relationship with both cluster recognition methods for fragments with 
Z = 3 — 8. In Fig.4, as an example, we plot the R21 = Y 2 /Yi values as a function of N of 
the emitted fragments obtained with the iso-MST and the MST cases. Here, "2" represents 
the neutron-rich reaction system 124 Sn + 124 Sn, and "1" represents 112 Sn + 112 Sn. The 
upper panel shows the results obtained with the iso-MST case and the lower panel shows 



16} , we find that the isoscaling 



results obtained with the MST case. As in previous study 
relationship exists in non-equilibrium transport model. The best fit isoscaling parameters a 
are listed on the figure and in Tabled. The isoscaling parameter a obtained with the MST is 
0.19, and the a value obtained with the iso-MST case increases to 0.22 which is still less than 
the experimental value of a = 0.36 ± 0.04 3l|. Sequential decays may modify the alpha 



values. However, results from current sequential decay calculations are model dependent 



32 



331 ] . Since the differences in a values are small, it will be difficult to determine the exact 



effects until more realistic sequential model calculations become available. 
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FIG. 4. (Color online) R21 values obtained from the ratios of the primary isotopic distributions for 
124 <Sn + 124 Sn collisions divided by those for 112 Sn + 112 Sn collisions, for 7^ = 0.5 at b=2fm. The 
upper panel shows the results obtained with the iso-MST case for Z=3 to 8 and the lower panel 
shows the results obtained with the MST case for Z=3 to 8. The best fit isoscaling parameters are 
listed in the panels. 
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Finally, we examine the effects of isospin cluster recognition on the system equilibrium 
measured by the production of light particles. To characterize the degree of equilibrium 
reached in the collision system, the observable of vartl = erf 



perimental and theoretical studies 



18 



19 



T lans Ing is often used in ex- 



35|. Here, the of mn and of are the variance of 



transverse and longitude rapidity distribution of fragment yields. When the reaction system 
reaches equilibrium, the value of crf rans /af ong is close to 1. In Fig. 5, we plot the calculated 
results of longitudinal (left panel) and transverse (right panel) rapidity distribution for the 
Z-weighted yield (Z=l-6) for 12A Sn + 124 Sn at b=2fm. The much wider longitudinal rapid- 
ity distributions suggest that the systems are far from equilibrium. More interestingly, the 
peaks in the longitudinal rapidity distribution (left pane) for Z=l-6 particles obtained in 
the MST case diminishes in the iso-MST case. It suggests that a higher equilibrium degree 
is predicted in the iso-MST case. With 7, = 0.5, the value of vartl = crf rans /af ong obtained 
with the iso-MST case is 0.62, which is slightly larger than 0.58 for the MST case. Thus 
in addition to the isospin sensitive observables, the equilibrium or stopping power of the 
system also depends on the detailed description of cluster formation implemented in the 
transport models as well as on the mean field and the in-medium NN cross section. 
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FIG. 5. (Color online) The calculated results of longitudinal (left panel) and transverse (right 
panel) rapidity distribution for the Z-weighted yield (Z=l-6) for 124 Sn + 124 Sn at b=2fm. 



In summary, we introduce a phenomenological isospin dependence in the description of 
cluster formation in transport models by adopting different Ro values for pp, nn and np, 
Rq P — 3/m and Rq 71 = Rq P = 6fm. Our results using the isospin dependent minimum 
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spanning tree method show suppression of Z=l particles and enhancement of fragments, 
especially for heavier fragments with Z > 12. Furthermore, we find enhanced production 
of neutron-rich isotopes at mid-rapidity. Consequently, isospin sensitive observables, such 
as the double ratios, DR(t/ 3 He), and isoscaling parameter a increase to larger values. The 
widths of the longitudinal and transverse rapidity distributions of Z = 1 — 6 particles also 
change. In all the observables that we examine, the effects introduced by the iso-MST 
algorithm are relatively small but in the direction of better agreement with data. However, 
we have not included sequential decays which may modify the magnitude of the effects. 
Nonetheless, the isospin dependence of the cluster recognition can be easily implemented 
and should be included in nuclear transport models. 
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